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It was recently shown by Motes, Gilchrist, Dowling &; Rohde [PRL 113, 120501 (2014)] that a time- 
bin encoded fiber-loop architecture can implement an arbitrary passive linear optics transformation. 
This was shown in the case of an ideal scheme whereby the architecture has no sources of error. 
In any realistic implementation, however, physical errors are present, which corrupt the output of 
the transformation. We investigate the dominant sources of error in this architecture — loss and 
mode-mismatch — and consider how it affects the BosonSampling protocol, a key application 
for passive linear optics. For our loss analysis we consider two major components that contribute 
to loss — fiber and switches — and calculate how this affects the success probability and fidelity 
of the device. Interestingly, we find that errors due to loss are not uniform (unique to time-bin 
encoding), which asymmetrically biases the implemented unitary. Thus, loss necessarily limits the 
class of unitaries that may be implemented, and therefore future implementations must prioritise 
minimising loss rates if arbitrary unitaries are to be implemented. Our formalism for mode-mismatch 
is generlized to account for various phenomenon that may cause mode-mismatch, but we focus on 
two — errors in fiber-loop lengths, and time-jitter of the photon source. These results provide a 
guideline for how well future experimental implementations might perform in light of these error 
mechanisms. 


I. INTRODUCTION 

In 2011 a simple model of quantum simulation us¬ 
ing linear optics - BosonSampling - was introduced 
by Aaronson & Arkhipov (AA) [ ]. This protocol is not 
believed to be capable of universal quantum computing 
as was shown possible in the seminal result by Knill, 
Laflamme & Milburn (KLM) [2] in 2001. Rather, Boson¬ 
Sampling can implement a subset of KLM linear optics 
quantum computing (LOQC) that does not require feed¬ 
forward, quantum memory, and dynamic control, only 
requiring single-photon Fock state preparation, passive 
linear optics, and photo-detection. BosonSampling has 
received much interest owing to its relative experimental 
simplicity as compared to universal quantum comput¬ 
ing, but is nonetheless believed to implement a classically 
hard problem. Several elementary experimental demon¬ 
strations of BosonSampling have recently been per¬ 
formed [3-7]. A major interest of late is to consider what 
other quantum states of light yield computationally in¬ 
teresting sampling problems that, like BosonSampling, 
yield a classically hard problem [8-T ]. 

BosonSampling is reminiscent of the Deutsch-Jozsa 
algorithm [12] of 1992. This algorithm demonstrated an 
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exponential separation between the complexity of clas¬ 
sical and quantum algorithms, but it solved a problem 
of no practical interest. Later, however, Shor described 
an efficient quantum factoring algorithm [13], a problem 
that is believed to be classically hard, generating much 
interest in the prospect of scalable quantum computing. 
Recently, and almost simultaneously, two separate works 
showed practical applications for BosonSampling. It 
was shown by Motes, Olson, Rabeux, Dowling, Olson & 
Rohde (MORDOR) that it may used for sub-shot-noise 
limited metrology [ ], and by Huh et al. that it can be 

used to generate molecular vibrational spectra [ ]. These 

applications are analogous too the Shor’s algorithm of 
the BosonSampling problem — a first glimpse into the 
potential of simple quantum sampling problems. For a 
more detailed introduction to BosonSampling see Ref. 
[ .6]. Although the fiber-loop scheme was initially pre¬ 
sented for the purposes of BosonSampling, Rohde re¬ 
cently demonstrated that with minor modifications the 
scheme can be made universal for quantum computing 
[ 7]. Here, however, we will focus on the application of 
this scheme to BosonSampling, or purely passive linear 
optics applications more generally. 

Although BosonSampling is much easier to imple¬ 
ment than universal LOQC, it remains experimentally 
challenging. The main challenges are preparing the input 
state and implementing the required unitary. It is possi¬ 
ble to prepare the desired Fock state input using sponta¬ 
neous parametric down-conversion (SPDC) single-photon 
sources [11, 18]. Loss in the unitary transformation is a 
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problem but there is evidence that even lossy systems or 
systems with mode-mismatch are still likely hard to sim¬ 
ulate given that the errors are sufficiently small [19, 20]. 
Constructing the required linear optics interferometer is 
challenging, as a Boson Sampling device might require 
thousands of optical elements, which must all be simul¬ 
taneously aligned. Two demonstrated ways to overcome 
the alignment problem are to use the time-bin encoded 
scheme by Motes, Gilchrist, Dowling & Rohde (MGDR) 
[ 51] or time-dependent dispersion techniques as presented 
by Pant & Englund [22]. Both methods do away with 
the hundreds or perhaps thousands of optical elements, 
requiring only a single pulsed photon-source and a sin¬ 
gle time-resolved photo-detector. An attractive feature 
of the former architecture is that there is only a single 
point of interference, and may therefore be much easier 
to align than conventional approaches. Additionally, the 
experimental complexity of these schemes are fixed, irre¬ 
spective of the size of the desired interferometer. In this 
manuscript we focus on analysing errors in the MGDR 
protocol to establish how well such a protocol might be¬ 
have in the presence of the two dominant sources of error 
— loss and mode-mismatch. 

MGDR showed that, using this architecture, arbitrary 
linear optics transformations can be implemented on a 
pulse-train of photons. However, this work assumes that 
the architecture has no sources of error. When errors are 
present the scheme no longer implements an arbitrary 
unitary transformation, but is constrained by the error 
model. In this work we analyse in detail various sources of 
error in the MGDR protocol. We begin by reviewing the 
MGDR architecture and then analyse the effects of lossy 
elements in the architecture followed by mode-mismatch 
caused by imperfect fiber-loop lengths and time-jitter in 
the source. These analyses accommodate the dominant 
challenges facing future experimental implementation. 


II. FIBER-LOOP ARCHITECTURE 


> rriT 



FIG. 1: The complete fiber-loop architecture fed by a 
pulse-train of photonic modes, each separated in time by 
r. The squares represent optical switches. A single length 
r inner fiber-loop is embedded inside a length > mr outer 
fiber-loop. The outer loop allows an arbitrary number 
of inner loops to be applied consecutively. When m — 1 
inner loops are implemented this architecture realises an 
arbitrary unitary transformation on m modes given that 
no loss is present. 


loop (after it traverses the inner loop once), such that the 
implemented unitary is bounded as an m x m matrix, 
where m is the length of the pulse-train. This can be 
obtained with, 

Ubs (1) - Ubs (m + 1) = (5 J ) > (!) 


where Ubs ( 2 ) * s th e unitary associated with the central 
beamsplitter at time i. The pulse-train then evolves in the 
inner loop. Each mode takes time r to traverse the inner 
loop so that it will interfere with the next time-bin at the 
central beamsplitter. Between each pulse a dynamically 
controlled beamsplitter, Ubs(^) of the form, 


UBs(t) 


( Ui'i (t) Ui t2 (t) \ 

V U2,l(t) U2, 2 (t) J ’ 


( 2 ) 


It was shown by MGDR that an arbitrary unitary 
can be implemented using time-bin encoding in a loop- 
based architecture [2'. ]. This is useful for the BosonSam- 
pling model because it significantly reduces the number 
of required optical elements. In fact, the experimental 
requirements to implement the architecture are fixed, ir¬ 
respective of the dimension of the desired unitary. Thus, 
the scheme is highly scalable, and uses far fewer physi¬ 
cal resources than schemes based on bulk-optics or inte¬ 
grated waveguides. In this architecture, shown in Fig. 1, 
a pulse-train of photonic modes consisting of, in general, 
Fock states and vacuum, are each separated by time r 
and sent into an embedded fiber-loop. The ith time-bin 
corresponds to the ith mode in a conventional spatially- 
encoded scheme. 

The boundary conditions of the protocol are that the 
first time-bin is coupled completely into the inner loop 
and the last time-bin coupled completely out of the inner 


where Ubs is an arbitrary, time-dependent SU(2) opera¬ 
tion, is applied. Here iqj is the amplitude of input mode 
i reaching output mode j. i = 1 (i = 2) represents the 
mode entering from the source (inner loop), and j = 1 
(j = 2) represents the mode exiting the loop (entering 
the loop). When a mode enters the loop it progresses to 
the next time-bin. 

After the entire pulse-train exits the inner loop the 
unitary map V is implemented, 

( 0 i > j + 1 

Vij = l U hl {i) i = j + 1 , 

{ ui )2 {i)u 2 ,i{j + 1) nu u 2 , 2 i k ) i<j + 1 

(3) 

where i G {1, m} and j G {1, m} represent input and out¬ 
put modes respectively. This map may easily be seen 
by writing out the spatial representation of the imple¬ 
mented unitary map and carefully following how each in¬ 
put mode traverses through to each output mode as done 
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in the original MGDR work. Note that we have employed 
a slightly different, but equivalent, indexing convention 
to the original MGDR proposal. 

The inner loop alone cannot implement an arbitrary 
unitary transformation so additional loops are required. 
The outer loop allows for an arbitrary number of appli¬ 
cations of the inner loop to be implemented. The net 
unitary U after L consecutive inner loops becomes, 

L 

U = l[V(l), (4) 

1=1 


where l denotes the Zth iteration of the inner loop. The 
pulse-train will traverse the inner loop L = m — 1 times 
and the outer loop m — 2 times for an arbitrary unitary 
transformation to be implemented. The outer loop must 
have round trip time > mr so that the pulse-train does 
not interfere with itself for a particular instance of the 
inner loop V(l). The pulse-train is coupled in and out 
of the outer loop via on/off switches. Once the desired 
transformation is performed, the pulse-train exits both 
loops and is measured via time-resolved photo-detection, 
where the time-resolution of the detector must be greater 
than r. The j th time-bin at the output corresponds to 
the jth spatial mode in the standard BosonSampling 
model. 

The unitary derived above assumes lossless compo¬ 
nents and perfect mode-matching at the central beam¬ 
splitter. In any realistic implementation this will not be 
the case, which we consider next. 


III. LOSS ERRORS 


In an implementation of a passive linear optics net¬ 
work, whereby the loss between each input/output pair of 
modes is uniform, loss simply amounts to a reduced suc¬ 
cess probability upon post-selecting on detecting all pho¬ 
tons. In the fiber-loop architecture, the different paths 
traverse the inner loop a different number of times lead¬ 
ing to non-uniform loss. This biases the unitary trans¬ 
formation resulting in a unitary that is not the desired 
one, even after post-selecting upon measuring all pho¬ 
tons. That is, the effects of loss cannot be simply fac¬ 
tored out of the unitary. In some architectures, asymmet¬ 
ric losses may be compensated for by artificially adding 
losses that rebalance the circuit, at the expense of over¬ 
all success probability. In the fiber-loop architecture this 
turns out not to be the case. 

In Sec. Ill A we introduce the metrics that we will use 
to analyse loss. In Sec. Ill B we determine the effect of loss 
due to the lossy switch and lossy fiber in the inner loop. 
Then in Sec. Ill C we analyse the net loss combining the 
inner loop losses with the outer loop losses. We denote 
quantities here that have loss with a prime. 


A. Loss Metrics 


1 . Similarity 


An interesting question is how small does loss need to 
be such that a particular unitary transformation is im¬ 
plemented with a particular error bar. The answer to this 
question is highly dependent on which unitary we wish 
to implement — some unitaries will suffer more asym¬ 
metric bias than others, depending on the switching se¬ 
quence that is required to implement them. Thus, the 
first question to ask is which unitary to consider. In the 
work of MGDR, a so-called ‘uniform’ unitary was con¬ 
sidered. This is a unitary where the amplitude (but not 
necessarily phases) of each element of the unitary are 
equal. That is, the magnitude of the amplitude between 
each input/output pair of modes is the same. This class 
of unitaries was considered as an example of ‘non-trivial’ 
matrices, which uniformly mix every input mode with 
every output mode. However, it is still an open question 
as to exactly what classes of unitaries yield hard sam¬ 
pling problems in the context of BosonSampling. We 
will here consider the same setting. We will explore this 
by using the similarity metric, <S, which compares the 
implemented map with the uniform map, 

(v;:,v o,, 2 - wa/ 

(E™=ii^i 2 )-(Er J =iiw iJ i 2 ) 

i fe-. Ml , 5) 

where IVij is an m x m uniform unitary given by 
\Wij\ 2 = 1/m. S is maximised by performing a Monte- 
Carlo simulation over different beamsplitter ratios so as 
to find the optimal switching sequence to make the map 
as uniform as possible. 


S = max 

Ubs (t) V t 


= max 
UBs(t) V t 


2. Post-selection Probability 

Another interesting question is how the probability of 
post-selecting upon all n photons is affected by loss, i.e 
the total success probability of the device. This is of espe¬ 
cial importance experimentally, as it directly translates 
to count rates. The post-selection probability of detecting 
all n photons at the output is, 

m (m \ ki 

*s»n x>«i 2 - ( 6 ) 

i=i \j=i ) 

where {k} is an integer string of length m that repre¬ 
sents a known input configuration of photons and ki is 
the number of photons in mode i. This equation is intu¬ 
itively derived as follows. For a single photon the prob¬ 
ability of entering mode i and exiting mode j is \Uij\ 2 . 
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Then the total probability that the ith photon exits the 
architecture is the sum of this over all j possible output 
ports, i.e. |^,j| 2 - Thus the probability of detecting 

all n photons at the output beginning in a particular con¬ 
figuration {k} is the product of this probability over all 
modes i where ki 7 ^ 0, as per Eq. 6 . This generalisation, 
by allowing arbitrary strings {&}, allows for implementa¬ 
tions such as randomised BosonSampling as described 
by Lund et al [11]. 

With losses present, U is in general no longer unitary. 
Rather, it is a mapping of input-to-output amplitudes, 
and will not be normalised. When there is no loss in the 
architecture Vs = 1? and with loss strictly Vs < 1, drop¬ 
ping exponentially with the number of photons. Imple¬ 
menting the required m — 1 loops will have exponentially 
worse loss than a single loop. 


B. Inner Loop Loss 



FIG. 2: A lossy inner fiber-loop fed by a pulse-train of 
photonic modes, each separated in time by r. We model 
the loss of the loop with a beamsplitter of reflectivity 
rjf and the loss of the switch with an efficiency r] s . Each 
mode experiences different amounts of loss, i.e. the first 
mode traverses the loop up to m times, the second up to 
m — 1 times, ..., and the 771 th mode at most once. 


C. Outer Loop Loss 

In the full fiber-loop architecture L inner loops are im¬ 
plemented via L — 1 round-trips of the outer loop, be¬ 
fore being coupled out to the detector. This architecture 
can implement an arbitrary unitary transformation when 
L = m — 1 if there are no errors present. The outer loop 
and outer switches cause a uniform loss on the entire 
pulse-train, since every path through the interferometer 
passes through these elements the same number of times. 
Hence, these factor out of U. The full lossy transforma¬ 
tion that occurs is then, 

u 1 = v ( l_i) n ^'(o. (8) 

1=1 

where L = m — 1 if an arbitrary transformation is de¬ 
sired, and V' is given by Eq. 7. The occurs 

because the pulse-train traverses an mr length of fiber 
in the outer loop L — 1 times (i.e rjf can be regarded 
as the efficiency per unit of fiber of length r), and the 
^ 2 (L-i) occurs because the pulse-train passes through 
the two outer switches L — 1 times. Fig. 3 shows the en¬ 
tire architecture with these loss errors. For an example 
of loop bias due to loss see App. A. Extending from this 
loop bias example we generalize the loss matrix denoted 
as £, which represents the accumulation of losses in the 
fiber-loop architecture, as a function of the number of 
loops L for an arbitrarily sized m x m transformation, 

C i j(L)=rtr) L+ i-\ ( 9 ) 

again where 77 = rjfrj s . Now the lossy map U' may be 
written as an element wise product of the ideal unitary 
U and the loss matrix £, 

U' = UojC. (10) 


We will model loss inside of the inner fiber-loop with 
a beamsplitter of reflectivity rjf and loss in the switch 
as 77s as shown in Fig. 2. When rjf = rj s = 1 the device 
has perfect efficiency. Before and after the inner loop the 
loss experienced by each mode in the fiber is negligible, 
since it may be arbitrarily short. Taking these losses into 
account, the implemented map of Eq. 3 becomes, 

0 i > j T 1 

Ui,i (i) i = j +1 

+ 1)- ’ ^ ' 

nUi^w i<j+ 1 

for a given loop, where 77 = r]fT] s . Note that this mapping 
is no longer a unitary matrix when 77/ < 1 or 77s < 1. This 
uneven distribution of losses in the input-to-output map¬ 
ping causes a skew in the matrix which prevents it from 
implementing the desired unitary transformation, even 
after post-selection. 


K> = ^ 1 


Elements of C that have no losses in them due to input 
modes not reaching output modes when L < m — 1 will 
be accounted for appropriately when C is multiplied by 
U by making the cooresponding matrix element in U' go 
to zero. 

In Fig. 4a we show how the optimised similarity 
with the uniform distribution varies with rjf and m for 
L = m — 1 inner loops, one photon in all m modes, and 
77 s = 1. With low loss rates (rjf « 1) the implemented 
unitary remains highly uniform. However, with several 
loops the success probability of detecting all n photons at 
the output decays exponentially as shown is Fig. 4b. For 
these plots the randomly generated U r that maximises S 
for each rjf and m is used to calculate the corresponding 

7>s. 

Now we consider how S and Vs are affected in Fig. 5 
with both the fiber loss and switch loss. We show this for 
the case of m = 3 and one photon per input mode, which 
is in the regime of present-day demonstrations. 
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FIG. 3: The full architecture which implements the lossy 
transformation U'. Each mode experiences 77 / m loss per 
outer loop since they each take time mr to traverse the 
outer loop. For an arbitrary unitary to be implemented 
in the ideal case the photons will traverse the outer loop 
L —1 times. This yields a net fiber loss from the outer loop 
of that can be factored out of J7', since it affects 

all paths equally. The net switch loss from the outer two 
switches is r] s 2 ^ L ~ 1 ^ and can also be factored out of U'. 
The losses within the inner loop, on the other hand, af¬ 
fect different paths differently, and in general cannot be 
factored out. 


IV. MODE-MATCHING ERRORS 

In any interferometric experiment it is inevitable that 
mode-mismatch will occur and is thus an essential source 
of error that we will consider in this section. There are 
many factors that may contribute to mode-mismatch in 
this architecture, such as incorrect fiber lengths, time- 
jitter in the sources, beamsplitter misalignment, and dis¬ 
persion of the wave-packets. In this section we will focus 
on two major sources of mode-mismatch: incorrect fiber 
lengths and source time-jitter. The former results in re¬ 
duced Hong-Ou-Mandel visibility at the central beam¬ 
splitter, owing to mismatched arrival times of photons. 
The latter effectively results in randomisation of the 
preparation times of the photons. 

We consider how mode-mismatch affects our protocol 
by calculating the fidelity, J 7 , between the ideal output 
state I'lpi) that one expects theoretically with no errors 
present, and the actual experimentally obtained output 
state | ipa). Imperfect fiber lengths and time-jitter both 
cause temporal shifts in the centre of the wave-packet, 
which will affect the output by both introducing uncer¬ 
tainty into the timing of the bins reaching the detec¬ 
tor, and undermining the Hong-Ou-Mandel visibility at 
the central beamsplitter. To calculate T then we need 
to calculate the temporal overlap between | xpi) and \^ a )- 
Therefore, we need to consider the temporal structure of 
the photons. 

We will model the temporal structure of photons using 
the formalism of Rohde et al. [23] . We only consider the 
inner loop in this analysis because there is no interfer¬ 
ence at any point in the outer loop. We obtain lower and 
upper bounds on T by performing a Monte-Carlo search 
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(b) 


FIG. 4: (a) Similarity S versus mode/photon number 
m = n and loop efficiency rjf for m — 1 loops. The map 
remains similar to the uniform unitary for low loss rates, 
implying that non-trivial unitary transformations may 
be implemented, m — 1 loops are considered because this 
is the number of loops required to implement an ar¬ 
bitrary unitary transformation in the lossless case, (b) 
Post-selection probability Vs versus mode/photon num¬ 
ber m = n and loop efficiency rjf for m — 1 loops. These 
two plots are related in that each point in Vs was cal¬ 
culated from the switching sequence U' corresponding to 
that which maximises <S. In both (a) and (b) the data 
was averaged over 1750 Monte-Carlo iterations and we 
let T] s = 1 , i.e the switches are ideal but the fibers are 
not. 


over different randomly generated unitaries V. We could 
also instead consider V' in this formalism to also jointly 
include losses. But we will treat losses separately from 
mode-mismatch for simplicity. 


A. Temporal Structure of Photons 

The temporal structure of a photon can be represented 
using a mode operator, 

/ oo 

— t — A)a^(x)dx, ( 11 ) 

-oo 
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FIG. 5: (a) Similarity <S, and (b) post-selection proba¬ 
bility Vs versus loop efficiency rjf and switch efficiency 
T] s with m = 3 modes, one photon per input mode, and 
m — 1 loops. These two plots are again related in that Vs 
is calculated from the switching sequence that maximises 
<S. This data was averaged over 1750 iterations. 


where ip(t — A) is the temporal density function centered 
at time £, A is a shift of the temporal centre of the pho¬ 
ton, and a) ( t ) is the time-dependent photon creation op¬ 
erator. This operator jV(t, A) acts on the vacuum |0) to 
create a photon with normalised Gaussian spectral den¬ 
sity function, 

^(x) = -jL=e~^, (12) 

V c V 7r 

where c/y/2 is the standard deviation. We assume that 
r A, in which case t denotes a time-bin, and A de¬ 
notes a small mismatch within the respective time-bin, 
not large enough to cause a photon to ‘jump’ from one 
time-bin to the next. Thus, both t and A represent shifts 
in the centre of the photon’s wavepacket, but the former 
is of the order of the time-bin separation, while the latter 
is of much smaller order than the time-bin separation. 

B. Our Formalism 

To analyse mode-mismatch we will consider three re¬ 
gions of the architecture we label as A, B , and C as shown 


in Fig. 6. Region A corresponds to the modes that are 
input into the architecture from the source, region B cor¬ 
responds to pulses inside the inner loop, and region C 
corresponds to pulses that exit the dynamic beamsplit- 


ter towards the detector. We introduce mode operators 
associated with each of these distinct regions — A^(t, A), 
B\t, A), and A) — each of the form of Eq. 11. 


T 


Region A 

m ••• 2 1 

1 -A A A 

.Reg,„B | E<!gi " C 

w > 


T 

UbsW) 



FIG. 6: The three regions we consider in the mode- 
mismatch formalism. Region A corresponds to the modes 
coming from the source, region B to the modes inside the 
inner loop, and region C to the modes exiting the loop. 


Since every pulse begins in region A the input state is 
a tensor product of pure states of the form, 


l*in> =' 






(13) 


where the tensor product is taken over all m modes, {k} 
is a known string representing the input photon-number 
configuration, and ki is number of photons in the ith 
input mode. 

Next, the input state is transformed by the dynamic 
beamsplitter, which takes the mode-operators from re¬ 
gion A into superpositions of regions B and C, 

UBs{t)A ] (t,A)ul s (t) -► u 1:2 {t)B^(t + 1, A) 

+ tii f i(t)Ct(t,A), (14) 

and pulses from region B to superpositions of regions B 
and C, 

U BS (t)B\t,A)Ul s (t) u 2>2 (i)B t (i + 1,A) 

+ A), (15) 

where we have used Eq. 2 for the elements of the dynamic 
beamsplitter at time t. Uusit) only acts on photons arriv¬ 
ing at the beamsplitter at time t =b A since r A. When 
a photon enters the loop t t + 1 as it advances to the 
next time-bin and will interfere with the next temporal 
mode. After this evolution, the entire pulse-train is cou¬ 
pled out of the loop such that the entire output state is 
a superposition of all possible output configurations. 

Now we model the state of the pulse train after t beam¬ 
splitters have been implemented, 
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n 

ji' = 1 

n ubsw 

_i'=1 


l*in> 



where the integer values of t denote the distinct time- 
bins. We note that there are ra + 1 total beam-splitters in 
a single implementation of the inner loop since there are 
771 — 1 beamsplitters to interfere the modes and another 
two beamsplitters to account for the initial and final 
boundary conditions of the MGDR protocol. Given how 
we modelled how the mode operators are transformed by 
C/bs W in Eqs. 14 and 15 the tth beam-splitter acts on the 
mode operators only in modal position t. Since a pulse 
coming out of the inner loop exits at beam-splitter t its 
modal position is m = t — 1 which accounts for there 
being m -hi beam-splitters and m modes. 

In general the final evaluated form of |T out ) may be ex¬ 
pressed as a superposition of all possible output photon- 
number configurations S', and their associated temporal 
configurations T(S), 


l*out> = £ £ 7S,TRC t (t Si , A r(s .)) 
S T(S) L i-i 


| 0 )® m , ( 17 ) 


where 7 s,t is the probability amplitude associated with 
photon time-bin configuration S and temporal shift con¬ 
figuration T(S), tSi denotes the time-bin of the ith pho¬ 
ton, T(S) denotes a configuration of temporal shifts asso¬ 
ciated with the configuration S, and A T(Si) is the tempo¬ 
ral shift of the ith photon associated with configurations 


A\u,Ai) k 


n 


y=l 


10 ) 


0 m 


(16) 


S and T. This is the most general representation of a 
configuration of photons across time-bins with associated 
shifts. The probability of measuring a particular config¬ 
uration is | 7 s 5 t| 2 , and to evaluate these probabilities we 
must fully characterise spectrum of time-bin and tempo¬ 
ral shift configurations, S' and T. Finding analytic forms 
for these expressions is largely prohibitive, and we calcu¬ 
late the 7 s,t via brute-force simulation of the evolution 
of the mode-operators through the network as described 
earlier. 


C. Fidelity Metric 

We analyse the results of this section by calculating the 
fidelity T between the ideal output state and the actual 
output state, given by, 

F=m\* a .)\ 2 , as) 

where |\Eq) is the ideal output state with no mode- 
mismatch (A 0) and |T a ) is the actual output state 
obtained with mode-mismatch. |T a ) reduces to |\tq) in 
the limit of no errors yielding T = 1. Calculating this 
overlap but letting |\fq) have general temporal mode mis¬ 
match until the end of the calculation we obtain, 


T = 


S' T'(S') 






S T(S ) 


i =1 


|0> 


0m 


E E 

S',S, T'(S'),T(S) L 


<*i| l^a) 

m m 

7s , ,T'7s,T(O| 0?n ]^[ C(ts> , A t >(s',)) A T ( S .))|O) 0m 


i '=1 


i=1 


(19) 


To simplify this expression further we use the formal¬ 
ism of second quantisation [2 ], which describes how 
the indistinguishability of particles in quantum me¬ 
chanics undergo symmetrisation. Here we use the ex¬ 
change symmetry of the bosonic Fock states, which ac¬ 
counts for how each temporal photon annihilation oper¬ 


ator C(ts', , At'(s'.,)) overlaps with each temporal pho¬ 
ton creation operator (T (£ 5 ., A T (£.)). Using bosonic ex¬ 
change symmetry we sum over all ml permutations of 
07 1 C(^,,A s 70™ 1 Ct(* Si ,As7 Then Eq. 19 be- 
comes, 




























T = 


E E 

S',S T'(S'),T(S) 


lS',T'lS,T ^ 


, At'cs^))^ (^, A T ( 5 .))|0) 


i=l 


( 20 ) 


where <j are the permutations over m elements. 

Finally, to calculate T we must find the wave packet 
simplification for (0| C(t s > a ,A T , (5 ^ ))C f (t^, A T(S .) |0), 

2' 7/ 


which we perform in App. B. Using this result we ob¬ 
tain, 


T = 


E E 

S',S T'(S'),T(S ) 


lS',T'lS,T 


E 


ri ex p( - 


( A r'(Si„) - a t(5 ; 


)f 


i=l 


4c 2 


( 21 ) 


Letting the ideal state |\Eq) have no temporal shifts, 
A 0, this reduces to, 


T = 


E E 7 S',T' ■ 75,t * ra! • exp 

5',5 T'(S'),T(S ) 


mA? 


T(5<) 


4c 2 


( 22 ) 


This derivation assumes the width of all wave-packets 
remain the same, i.e the photons are identical up to a 
temporal displacement. The width of the wave-packets 
may broaden due to dispersion but under the relatively 
short lengths of fiber-loop required for small m the effect 
of dispersion may be neglected; however this formalism 
may be easily modified to include dispersion by creating 
an operator that broadens the wave-packet width c as a 
function of the length of the fiber the wave-packet has 
traversed. 

Next we consider two types of mode-mismatch: non¬ 
ideal lengths of the inner loop, and time-jitter at the in¬ 
put source. 


r + 5 



FIG. 7: The inner fiber-loop with an error S in its in¬ 
tended length. Every time a pulse traverses the inner 
loop it is shifted from its expected temporal position by S , 
thereby reducing the Hong-Ou-Mandel visibility at Ubs- 


D. Imperfect Inner Loop Length 

Here we analyse errors in the MGDR fiber-loop archi¬ 
tecture caused by a non-ideal length of inner fiber-loop as 
shown in Fig. 7. We let the length of the inner loop have 
some length r + S, where 6 is the error in the intended 
length r and may be positive or negative. Thus every 
photon that traverses the inner loop acquires a temporal 
shift of S from its expected centre. We ignore imperfect 
lengths of the outer loop because every mode will tra¬ 
verse the outer loop an equal number of times creating a 


global temporal shift with no impact on interference at 
the central beamsplitter. 

The input state is given by Eq. 13 where A^ = 0 V i. 
This models an ideal input state with no time-jitter or 
other errors in the source. To account for the unwanted 
time-delay S we introduce the time-delay operator T(5), 

T(6)B\t, A)T t (<5) = &{t, A + 6), (23) 

which acts only in region B - the region inside the in¬ 
ner loop. This adds a small temporal displacement, not 
enough to confuse time-bins. Thus it affects A but not 
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t. It has no effect on the mode-operators A and C. Using 
the boundary conditions shown in the MGDR protocol, 
the first photon is coupled completely into the loop so 
it picks up a time-delay of S. Afterwards the pulse-train 
interacts at the beamsplitter described in Eqs. 14 and 
15, where it is sent into a superposition of regions B and 


C. As the state evolves all amplitudes entering the inner 
loop (region B) will acquire a time-shift of S. After the 
last mode traverses the inner loop the state is coupled 
completely out as per the MGDR protocol. The output 
state is given by, 




n twiWi') 

_i' = 1 


m 1 


ki 


n twcwo 


U'=l 


|o> 


<S>m 


(24) 


where we have inserted the time-delay operator appropri¬ 
ately in Eq. 16. Fig. 8a shows how the fidelity T scales 
with m, n, and S and Fig. 9b shows the worst- and best- 
case fidelities, where we have searched over switching se¬ 
quences. 

E. Time-jitter from Input Source 

A major source of error in the time-bin architecture is 
time-jitter of the input source. Ideally each mode will be 
separated by time r but in reality non-ideal sources will 
randomly shift modes from their desired centre of time U 
in mode i. To model time-jitter we let the temporal shift 
of input mode z be a Gaussian random variable drawn 
from the normal distribution, 

m, - ] = A? xp • (25) 

centered in mode i at time U and with a standard devi¬ 
ation of a. The input state of Eq. 13 becomes, 


We assume that the shifts caused by time-jitter are 
much less than the time-bin separation r, such that the 
probability of time-bin confusion remains negligible, i.e. 
A f(ei) <r. Fig. 9a shows how the fidelity T scales with 
m, n, and a. Fig. 9b shows the worst- and best-case J 7 , 
searching over many switching sequences. 

V. CONCLUSION 

In this work we have analysed sources of error in the 
Motes, Gilchrist, Dowling & Rohde fiber-loop architec¬ 
ture for implementing BosonSampling. Specifically we 
have analysed loss and mode-mismatch. In the loss anal¬ 
ysis we examined how lossy fibers and switches affect the 
operation of the architecture in both the inner and outer 
loops. We found that loss causes an asymmetric bias in 




2.0 


(b) 

FIG. 8: (a) The average fidelity T between the ideal state 
|\Eq) and the actual experimental state |T a ) versus the 
error in the intended length of the inner loop 5. (b) The 
worst (bottom) and best (top) case fidelity T between 
the ideal state |\Eq) and the actual experimental state 
|T a ) versus the error in the intended length of the inner 
loop S and number of modes m. In (a) and (b) there are 
m modes with one photon per mode and the data was 
obtained over 250 implementations each with a unique 
randomly generated unitary. 


the desired unitary, unique to a temporally implemented 
unitary transformation. That is, even upon post-selection 
the operation of the device is erroneous. Additionally, 


|0}i. (26) 


IV’in) = 0 
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(b) 


FIG. 9: (a) The fidelity T between the ideal state |\fq) and 
the actual experimental state |T a ) with random time- 
jitter in the input source versus modes m and standard 
deviation a with no fiber length error S = 0. (b) The 
worst (bottom) and best (top) case fidelity T between the 
ideal state |\Eq) and the actual experimental state |T a ) 
with time-jitter. In (a) and (b) there is one photon per 
mode, the data was averaged over 250 implementations 
each with a unique randomly generated unitary, and the 
time-jitter was drawn from the normal distribution. 


like all linear optical architectures, our scheme has ex¬ 
ponential dependence on loss, thereby reducing the post¬ 
selection success probability of detecting all n photons. In 
the mode-mismatch analysis we analysed only the inner 
loop since no interference occurs in the outer loop. We 
examined two types of mode-mismatch including an im¬ 
perfect length of fiber in the inner loop, and time-jitter of 
the photon source. This analysis provides a guideline for 
future experimental implementations, to provide insight 
into how such a device might realistically behave in the 
presence of loss and mode-mismatch, the two dominant 
error mechanisms affecting this protocol. 
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traversal of the outer loop), we find the input-to-output 
loss matrix to be, 


( Vs^f 2 r] s 5 r]f 3 
V Vs 3 Vf Vs i Vf 2 



(A2) 


where we have ignored the losses due to the outer loop 
as it fields an overall normalisation factor that does not 
bias U. As we can see, for each iteration of the inner 
loop U accumulates more loss, with a decreasing overall 
success probability, but the amount of skew in the matrix 
remains the same. 


Appendix B: Wave-packet Simplifications 


Appendix A: Intuitive Example of 
Loop Bias due to Loss 

An example of how U becomes biased is explained here. 
Let’s consider two examples of a two-mode pulse-train — 
a single inner loop and two inner loops. 

1. One Loop 


In this section we derive the overlap of two photons 
with temporal creation operator Af(t, A) as per Eq. 11 
and temporal density function A). For the purpose of 
this work we assume that the temporal spacing r between 
each mode is much larger than the width of the wave 
packet c such that the overlap of our temporal photons 
in different time-bins is negligible, 

(0|i(£,A)it(r,A')|0) =0, (Bl) 


The first mode can can exit the first output mode by 
traversing the inner loop once. Here it picks up loss due 
to the middle switch twice rj s 2 , and loss due to the inner 
loop fiber once r)f, obtaining a net loss of r] s 2 r]f. The 
first mode can exit the second output mode by traversing 
the inner loop twice. In this case it obtains a net loss of 
r] s 3 r]f 2 . A similar analysis can be performed for the other 
combinations. Then, we can write the loss amplitudes 
corresponding to the input (rows) and output (columns) 
modes in matrix form as, 


( r ls 2 r]f Vs 3 Vf 2 

V Vs vs 2 vf 



(Al) 


where r] = r] s r]f and observe the bias accumulating in this 
input to output map. The net input-to-output mapping 
of amplitudes is given by taking the element-wise product 
of this loss matrix with the ideal unitary, £oL, thereby 
leaving us with a biased map. 


for t 7 ^ t'. 

For photons in the same time-bin and allowing for ar¬ 
bitrary temporal-shifts, the overlap is, 


(0|A(£, A')A\t, A)| 0 ) 

f poo 

(0| / if*(x f — t — A')a(x')di 

K J — GO 

J ip(x — t — A)a^ (x)dx\0) S j 

/ OO POO 

/ x/;*(x f — t — A')ip(x — t — A) 

-oo J —oo 

x (0|a(x')d^(x)|0) dx'dx 


/ oo 

ip*(x — A')ip(x — A )dx 

-oo 


-oo 
(A 7 —A ) 2 

= e \ 


(B2) 


2. Two Loops 

A similar analysis as above but following the paths for 
two consecutive applications of the inner loop (i.e one 


For ideal states where A' = A = 0, we notice that 
J 7 = 1, as expected when there is no mode-mismatch. We 
use these results for simplifying Eq. 21 in our analysis of 
mode-mismatch. 




